Levinson's Theorem for the Klein-Gordon Equation 

in Two Dimensions 



Shi-Hai Dong* 

Institute of High Energy Physics, P. O. Box 918(4), Beijing 100039, The People's Republic of China 

Xi-Wen Hou 

Institute of High Energy Physics, P. O. Box 918(4), Beijing 100039 
and Department of Physics, University of Three Gorges, Yichang 443000, The People's Republic of China 

Zhong-Qi Ma 

China Center for Advanced Science and Technology (World Laboratory), P. O. Box 8730, Beijing 100080 
and Institute of High Energy Physics, P. O. Box 918(4) Beijing 100039, The People's Republic of China 

Abstract 

The two-dimensional Levinson theorem for the Klein-Gordon equation with 
a cylindrically symmetric potential V(r) is established. It is shown that N m ir = 
77 ( n m ~ n m) = [àm(M) + fa] — [S m (—M) + fa] ,where N m denotes the difference 
between the number of bound states of the particle n+ and the ones of antiparticle 
n~ with a fixed angular momentum m, and the 5 m is named phase shifts. The 
constants fa and fa are introduced to symbol the critical cases where the half 
bound states occur at E = ±M. 
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I. Introduction 



The Levinson theorem[3] , an important theorem in scattering theory, established 
the relation between the total number of bound states and the phase shifts at zero 
momentum. During the past half century, the Levinson theorem has been proved by 
several authors with different methods, and generalized to different fields [4-11]. Rough 
speaking, there are three main methods used to prové the Levinson theorem. One [3] is 
based on the elaborate analysis of the Jost function first introduced by Jost. The second 
is relied on the Green function method [7]. The third method is used to demonstrate 
the Levinson theorem by the Sturm-Liouville theorem [8-10]. This simple, intuitive 
method is readily to be generalized and has been verified by the proofs of many physical 
problems [8-10,22-24]. Furthmore, some obstacles and ambiguities, which may occur 
in other two methods, disappear in the third method. However, it is found in the later 
proof that the Sturm-Liouville theorem can't be directly used to prové the Levinson 
theorem for the Klein-Gordon equation, but a modified method which is similar to the 
Sturm-Liouville theorem will be applied to prové the Levinson theorem. Consequently, 
such a generalization may be useful for the method of bosonization method which has 
been widely utilized in the literature[26]. 

The Klein-Gordon equation, which describes the motion of a relativistic scalar par- 
ticle, is a second-order differential equation with respect to both space and time. When 
there exists a potential as the fourth component of the vector field, the energy eigen- 
values are not necessarily real and the eigenfunctions satisfy the orthogonal relations 
with a weight factor [1-2] such that a parameter e which is not always real and positive 
appears in the normalized relation with a weight factor. As pointed out by Pauli and 
Snyder at al[l-2] ,after bose quantization, that those amplitudes with real and positive 
e describe partides, but those with real and negative e antiparticles. 

Recalling in the three-dimensional spaces, two main methods are used to set up 
the Levinson theorem for the Klein-Gordon equation. One is relied on some formulae 
which are vàlid for the cases without complex energies [7]. The other, which is similar 
to that of Sturm-Liouville theorem, is applied to arrive at the Levinson theorem for the 
Klein-Gordon equation[9]. This result is correct for the cases both without complex 
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energies and with complex energies. 

The reasons why we write this paper are that, on the one hand the Levinson theorem 
in two dimensions has been studied in experiment [19] as well as in theory [20-24] in 
virtue of the wide interest in lower-dimensional field theories and other modern physics 
[12-18], on the other hand the Levinson theorem for the Klein-Gordon equation in two 
dimensions has never been appeared in the literature. In our previous works[22-24], 
some surprised results are obtained from the nonrelativistic and relativistic particle as 
well as the non-local interactions in two dimensions. We attempt to set up the Levinson 
theorem for the Klein-Gordon equation in two dimensions. What new results will be 
appeared? 

This paper is organized as follows. In Sec. II, we review the properties of the 
Klein-Gordon equation, especially those related with the parameter e. In Sec. III, it 
is proved that the difference between the numbers of bound states of particle and the 
ones of antiparticle only relies on the changes of the logarithmic derivatives of the wave 
functions at E = ±M as the potential V(r) changes from zero to the given value. In 
Sec. IV, it is also turned out that these changes are closely connected with the phase 
shifts at E — ±M which then results in the establishment of the two-dimensional 
Levinson theorem for the Klein-Gordon equation. 



Throughout this paper the natural units h — c — 1 are employed. Consider a 
relativistic scalar particle satisfying the Klein-Gordon equation 



where the potential V(x) is the fourth component of a vector field and the M, E denote 
the mass and the energy of the particle, respectively. In order to simplify the discussion, 
we only research that the potential is static and cylindrical symmetric one 



II. the Klein-Gordon Equation 




(i) 



V(x) = V(r), 



(2) 



and its asymptotic behavior is written 



r|V(r)| — > when r 



(3a) 
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and 

V(r) = when r > r . (36) 

Equation (3a) is required to make the wave function single value at the origin, and 
(3b) is called the cutoff potential for the sake of the simplicity of discussion, i.e it is 
vanishing beyond a sufficiently large radius r Q . It is proved that, following the method 
[22-23], the results obtained in this paper will not change the essence of the proof if 
the potential vanishes faster than r~ 2 at infmity. 
Introduce a parameter À for the potential V(r) 

V(r,\) = \V(r), (4) 

which shows that the potential V(r, À) changes from zero to the given potential V(r) 

when À increases from zero to one 

Due to the symmetry of the potential, Let 

V>(x, A) = r- 1 / 2 R m (r, A)e ±im ^, m = 0, 1, 2, . . . , (5) 

where the radial wave equation R m (r, À) satisfies the radial equation 

d l30R + _ M 2) _ (2EV _ ^ _ tL·AIí j flm(r , A) = 0. (6 ) 

Denote by R m i(r, A) the solution to Eq.(6) for the energy E 1 

d 23^hR + _ M ^ _ (2ElV _ V) - !^Í) JU(r, A) = 0. (7) 

Multiplying Eq.(6) and Eq.(7) by R m i( r , A) and R m ( r , A), respectively, and calculating 
their difference, we have 

|- {i? m (r, A)id(r, A) - i? ml (r, A)^(r, A)} = —{E{—E)R* ml (r , A) • (£? 1 *+£?-2y)i2 rïl (r, A), 

(8) 

where the primes denote the derivative of the radial wave function with respect to 
the variable r. As we know, the energy eigenvalues are not necessarily real for some 
potential V(r) which origins from the Klein paradox. Integrating (8) over the whole 
space and noting that R m (r, A)i?^(r, X) — R ml (r : A)i?^ t (r, A) vanishes both at the origin 
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and at infinity for the physically admissible solutions with the different energies E and 
Ei, we get the weighted orthogonality relation of the radial wave function 

{El - E) / iCi(r, \){E*i +E- 2V)R m (r, X)dr = 0. (9) 
Jo 

As a matter of fact, we always able to obtain the real solutions for the real energies. 
However, it is easy to see from Eq. (9) that the normalized relation for the solutions 
with real energies are not always positive on account of the weight factor (E^ + E — 2V): 

/■oc í e E S(E 1 - E), \E\ > M, 

/ R ml (r,X)(E 1 + E-2V)R m (r,X)dr= (10) 

Jo { e E 5 ElE , \E\ < M. 

The parameter ee, which depends on the particular radial wave function R m (r,X), 
may be either positive, negative or vanishing. Normalized factors of the solutions 
can't change the sign of e. Generally speaking, if the solution R m (r, X) with a complex 
energy E is complex, then is also a solution with complex energy E* and a complex 
ee appears for a pair of the complex solutions. It is evident after bose quantization 
that those R m (r,X) with positive ee describes partides and those with negative €e 
antiparticles. In the case zero ee, the solution can be regarded as a pair of particle and 
antiparticle bound states. The Hamiltonian and charge operator can't be written as the 
diagonal forms for the solutions with complex energy ee, therefore they describe neither 
partides nor antiparticles. In this paper, we only count the number of bound states 
with the real positive and negative nonvanishing e E is named particle and antiparticle 
bound states, respectively. 

Since we are always able to arrive at the real solution for the real energy, we can 
now solve Eq.(6) in two regions and match two solutions at r . Actually, the solutions 
in the region [0,r ] with R m (0) = can be arrived at in principle. We only need one 
matching condition at ro for the logarithmic derivative of the radial wave function 

(ii) 

Only one solution is convergent at the origin because of the condition (3a). For 
example, for the free particle (À = 0), the solution to Eq. (6) at the region [0, ro] is 
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proportional to the Bessel function J m (x): 



R m (r,0) = { V 2 



7T h'T' 

" J m (kr), when |£| > M and fe = ^E 2 - M 2 



e- ím7r /y^J m (mr), when |£| < M and « = y/M 2 - E 2 , 

(12) 

The solution R m (r, 0) given in Eq. (12) is a real function. A constant factor on the 
radial wave function R m (r,0) is not important. 

In the region [r ,oo), we have V(r) = 0. For \E\ > M, there are two oscillatory 
solutions to Eq. (6). Their combination can always satisfy the matching condition 
(11), so that there is a continuous spectrum for \E\ > M. 



R m (r, X) = ]l — {cosr] m (k, \)J m (kr) - smr] m (k, X)N m (kr)} 



~ cos 



(*r-^-ï + W*,A)), whenr^oo. (13) 

where N m (kr) is the Neumann function. 

However, there is only one convergent solution in the region [r , oo) for \E\ < M 
the matching condition (11) is not always satisfied. 



R m (r, A) = e i{m+1)n/2 ^l^-H^(ÍKr) ~ e~ Kr , when r — ► oo. (14) 

where H$(x) is the Hankel function of the first kind. When the condition (11) is 
satisfied, a bound state appears at this energy. It means that there is a discrete 
spectrum for \E\ < M. 

As mentioned above, in the case with the real energy solutions, integrating the Eq. (6) 
in two regions [0,r ] and [r ,oo) ,respectively, and taking the limit E\ — * E, we will 
obtain the following equations in terms of the boundary condition that R m (0) = and 
iC(oo) = for \E\ < M 

dA m (E,X) _ d ( 1 dR m (r,X) \ 

OE ~ dE\R m (r,X) Or ) r=rQ _ (15o) 

= -R m (r ,X)~ 2 R m (r,X) 2 2[E-V(r)]dr<0. 
Jo 

and 

dB m (E) _d_ ( 1 dR m (r,X) \ 

dE ~ dE\R m (r,X)^ dr ) r=r<j+ (15fe) 

= R m (r ,X)- 2 I" R m {r,X) 2 2 E dr > 0. 

Jro 
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which demonstrates from Eq. (15) that A m (E, A) is no longer monotonic with respect 

to energy, but B m (E) is still monotonic with respect to energy if the energy doesn't 

change sign. 

From the matching condition (11) we have 

. (h u_ Jm(kro) A m (E, A) - kJj(kr )/J m (kr ) - l/(2r ) 

Vm[ ' } N m (kr ) ' A m (E,X)-kN^kr )/N m (kr )-l/(2r o y ^ 

Vm(k) = r} m (k, 1). (17) 

where the prime denotes the derivative of the Bessel function, the Neumann function, 
and later the Hankel function with respect to their argument. However, it is not true 
for \E\ < M because of no adjustable phase shift S m (E). Once the matching condition 
is satisfied, we will get the discrete bound states. 

The phase shift r] m (k, X) is determined from (16) up to a múltiple of n due to the 
period of the tangent function. In this paper, for the free particle (V(r) = 0), the 
definition of phase shift i] m (k, 0) is defined to be zero, i.e 

r] m (k,0) = 0, where A = 0, (18) 

which is same as our previous definition [8- 9, 22-24]. 

It is shown from Eq. (10) that scattering states — -\E\ > M — are normalized as 
the Dirac S function, and that the main contribution to the integration Eq. (10) comes 
from the radial wave functions in the region[r , oo) where there is no potential. For 
this reason we obtain 

E 



e E = ttVE 2 - M 2 ■ — , \E\ > M. (19) 

\E\ 

All the scattering states with positive energy (E > M) describe partides and those 
with negative energy (E < —M) describe antiparticles. It is easy to see that this 
conclusion is not true for the critical case E = ±M except for S waves where there is 
a half bound state at E — ±M. The situations which e# with E = ±M and m > 1 
may be positive, negative or vanishing are relied on the potential. 

III. The Number of Bound States 



In our previous works, the Levinson theorem for the nonrelativistic and relativistic 
partides are set up under the help of Sturm-Liouville theorem. For the Sturm-Liouville 
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problem, the fundamental trick is the definition of a phase angle which is monotonic 
with respect to the energy [25]. Although this method is very simple, intuitive and easy 
to be generalized, from the Eq. (6), it is the weight factor (E{ + E — 2V) that makes the 
Sturm-Liouville theorem not be used for the the Klein-Gordon equation. Nevertheless, 
a modified method is applied to prové the Levinson theorem for the Klein-Gordon 
equation. From the difference between the Eq.(15a) and Eq. (15b), we arrive at 

dB m {E) 8A m {E,X) _ 1 

~dE dÈ— = Bm{E) ~ Am{E ' A) " BM €E > (20) 

where here and hereafter the primes denote the derivative with respect to the energy. 
From Eq. (14), we get 



m{ ) ~ Hrn {íktq) ~2r - ' 



(-m + l/2)/r = p m when k x — >0 

[21) 

—k ~ — oo when k\ — > oo. 



The logarithmic derivative given in Eq. (21) does not depend on À. On the other hand, 
when À = we obtain from Eq. (12) 



A m (E,0)-. mJ - (mro) 



J m (inr ) 2r 



(m + l/2)/r when ki — >0 

(22) 

k ~ oo when k\ — > oo. 



It is evident from the Eqs. (21) and (22) that both B m (E) and A m (E, 0) are continuous 
curves with respect to energy which don't intersect each other; i.e. the matching 
condition (11) is not satisfied if \E\ < M and À = 0. No bound states appear when 
there is no potential. 

As À changes from the zero to the given potential, B m (E) don't change, but 
A m (E,\) changes continuously except the points where R m (r ) = and A m (E, X) 
tends to infinity. Generally speaking, A m (E, X) is continuous except those finite points 
and intersects with the curve B m (E) several times for \E\ < M. The bound state will 
appear only if the intersection happens. The points of the intersection determine the 
number of the bound states. It is shown from Eq. (20) that the relative slopes at the 
points of intersection decide whether the bound states describe particle or antiparticles. 

When the potential V(r) change with the À, the number of intersection points will 
change, too. This only origins from the following two sources. Firstly, the intersection 
points move inward or outward at E — ±M. Secondly, the curve A m (E, X) intersects 



with the curve B m (E) or departs from it through the tangency point. For the second 
case, a pair of particle and antiparticle bound state will be created or annihilated at the 
same time, but the difference of the number of the particle n+ and antiparticle bound 
state n~ don't change. That's to say, the change of the whole bound states N m which 
expresses that the difference of the particle bound state and the antiparticle state only 
depends on the intersection points moving in or out at E — ±M where the critical cases 
occur. Hence, we only discuss this case. There are four cases when A m (M, X) = B m (M) 
when A m (M, X) decreases across the value B rn (M) = (— m + l/2)/r at E — M 

(l)A' m (M,X)<B' m (M), (23a) 

(2) A£\M,\) = B%\M), (-l) n A^ +1 \M,X) < (-l) n B% +1 \M), (236) 

(3) 4?(M, A) = BW(M), (-1)V^ +1 )(M, A) > (-l) n B^ +1 \M), (23c) 

(4)A' m (M,X)>B>(M). (23d) 

Where here and hereafter n is positive integer. For the first two situation, a inter- 
action point moves inward from E > M to E < M, which results in the appearance 
of new particle bound state. However, for the last two cases, the interaction point 
moves outward from E < M to E > M, which causes the disappearance of an an- 
tiparticle bound state. The converse process occurs when A m (M, X) increases to cross 
the value B m {M), i.e. the number of bound states N m increases by one only if each 
time A m (M, X) decreases to cross the value B m (M) at E = M. Conversely, each time 
A m (M, X) increases across the value B m (M) at E = M, N m decreases by one. 

On the other hand, there are also four cases when A m (— M, X) = B rn (—M): 

(ï)A' m (-M,X)>B' m (-M), (24a) 

(2>)A%\-M, X) = B^(-M),A^ +1 \-M,X) > ^ +1) (-M), (246) 

(3')A$(—M, X) = B^\-M),A^ n+1 \M,X) < B^ n+1 \M), (24c) 

(4')A' m (-M,X)<B' m (-M). (24d) 

If A m (— M, X) decreases across the value B m (—M) as À increases, for the first two 
cases a interaction moves inward from the E < —M to E > —M point which describes 
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an antiparticle. But for the last two cases an interaction point moves outward from 
E > —M to E < —M which describes a particle. The number of bound states N m 
decreases by one only if each time A m (—M, X) increases across the value B(—M). The 
opposite process occurs when A m (—M, X) increases across the value B rn (—M). 

We denote by N m (±M) the difference between the number of times A(±M, X) de- 
creasing across the value B(±M) and the number of the times that A(±M, X) increasing 
across that value. Hence, we obtain 

N m = n+ -n~ = n m (+M) -n m {-M). (25) 

IV. The Phase Shifts 

As we know, the solutions in the region [r , oo) for the scattering states have been 
given by Eq. (13). The phase shift f] m (Q, X) is the limit of the phase shift r] m (k,X) 
as k tends to zero. Hence, what we are interested in is the phase shift f] m (k, X) at a 
sufficiently small momentum k, k <C l/r . For the small momentum we obtain from 
the matching condition (11) 



tan. rj m (k, A) ~ 

-n(kr ) 2m A m (0, X) - (m + l/2)/r 



2 2 ™m!(m - 1)! ^(o, A) - c 2 k 2 - p m ( 1 ^ 



(m- l)(2m- 1) 
Akro) 2 A m (0, X) - 3/(2r ) 



4 A m (0, X) - c*k* - Pi (1 + 2(kr y \og(kr )) 

A m (0,X)-c 2 k 2 -p (l-(kr f) 



when m > 2 

when m = 1 
when m = 0. 



21og(A;ro) ( _2 \ 

(26) 

In addition to the leading terms, we include in (26) some next leading terms, which is 
useful only for the critical case where the leading terms are canceled with each other. 
and 



dr] m (k, X) 
dA m (E,X) 



-8r cos 2 t] m (k,X) 



tt {2r A m (E, X)N m (kr ) - 2kr N^(kr ) - N m (kr )Y 



< 0, (27) 
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which shows that the phase shift is monotonic with respect to the logarithmic derivative 
A m (E, À) as À increases. 

It is shown from Eqs. (26) and (27) that they are not different from those of 
Schròdinger equation. Therefore, we may simply discuss this problem by the same 
method. Each time A m (±M,\) decreases across the value B m (±M) as the potential 
changes from the zero to the given potential, the phase shift S m (±M, À) increases by ir. 
Conversely, the phase shift 8 m (±M, À) decreases by ir if A m (±M,\) increases across 
the value B m (±M). 

As À increases from zero to one, i.e. the potential changes from the zero to the 
given value, we have 

6 m (±M) = 6 m (±M, 1) = n m (±M)n, (28) 

Thus, we draw a conclusion that the Levinson theorem for the Klein-Gordon equation 
if A(±M, 1) yé B(±M) 

N m 7r = 5 m (M)-5 m (-M). (29) 
We now discuss the critical cases 

A m (M, 1) = B m (M) and A m (-M, 1) = B m (-M), (30) 

where the potential changes from the zero to the given potential V(r). Similar to the 
discussion [22-24], the phase shift 5 m (±M, À) increases by n for m > 1 or an additional n 
for the P waves if A m (±M, À) decreases from near and larger than the value B m (M) to 
smaller than that value when the potential changes from the zero to the given potential. 
Conversely, 5 m (±M, À) doesn't decrease by n or an additional n if A m (±M, À) increases 
across the value B m (M) as the potential changes to the given potential V(r). On the 
other hand, the states for m — 0, 1 are called a half bound state which is defined as its 
wave function is finite but not square integrable. Furthermore, the half bound state 
is not a bound state. For M > 1 states in the critical situations, there is a bound 
state but its €e may be either positive, negative or vanishing, which depends on the 
different cases (23) and (24). We consider the state with zero e^asa pair of particle 
and antiparticle bound states. 
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Introduce two parameters (3\, to describe the appearance or disappearance of the 
bound states at the critical cases. fi\ = for the noncritical case A(M, 1) ^ B m (M), 
and (3 2 = for the case A m (—M, 1) ^ (-M). 

(1) If A(M, 1) = B m (M),p 1 = for the cases (23a) or (23c) with m > = -1 for 
the cases (23b) or (23d) with m > 1; and = — 1 for the case (23a) with m—1. 

(2) If A m (—M, 1) = B m (-M),f3x = for the cases (24a) or (24c) with m > = -1 
for the cases (24b) or (24d) with m > 1; and /3i = — 1 for the case (24a) with m—1, 
where À is substituted by one. Then, the Levinson theorem for the Klein-Gordon 
equation with the cylindrical symmetric potential V(r) satisfying the asymptotic be- 
havior(3) 

N m n = vr (n+ - N~) = [6 m (M) + ft] - [5 m (-M) + /3 2 ]. (31) 

According to the above discussion, it is easy to find, compared with the case in 
the three-dimensional spaces, that the phase shifts for the critical sates changes by an 
additional n not by n/2. This conclusion is same as the relativistic and nonrelativistic 
partides. 
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